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^ I Abstract 

We calculate connected and disconnected contributions to the flavour singlet 
\q , scalar density amplitude of the nucleon in a full QCD lattice simulation with rif = 2 

dynamical Wilson fermions at /3 = 5.6 on a 16 3 x 32 lattice. We find that both 
contributions are of similar size at the light quark mass. We arrive at the estimate 
00 . ovat = 18(5)MeV. Its smallness is directly related to the apparent decrease of u, d 

quark masses when unquenching QCD lattice simulations. The y parameter can be 
estimated from a semi-quenched analysis, in which there are no strange quarks in the 
sea, the result being y = 0.59(13). 

1 Introduction 

> . 

X ■ The pion-nucleon er-term is defined as the flavour singlet scalar density amplitude of the 
nucleon, multiplied by the light quark mass m udi 

a nN = m ud (N\uu + dd\N) , m ud = l/2(m u + m d ) . (1) 

Our motivation to study this quantity in a full QCD simulation is twofold. First of all u^n 
provides a direct measure of the explicit chiral symmetry breaking in QCD. A comparison 
of its experimental value with a (first principles) QCD calculation is therefore of great im- 
portance for the understanding of the chiral properties of the strong interaction. Secondly, 
the scalar density amplitude receives contributions from both connected and disconnected 
(vacuum polarization) diagrams. This is shown in fig.[l]. In a flavour singlet combination, 
such as cr n N, the latter are simply added. Therefore, one expects sizeable contributions 
from those processes, which most likely will depend on the details of the vacuum structure. 
Given this, the pion-nucleon cr-term provides an ideal opportunity to disclose the impact 
of sea quarks on nucleon properties. 

On the other hand the determination of the experimental value of a n N is by no means 
straightforward as it requires quite a bit of theoretical input. For instance it implies the use 
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is related to cr^N by 



Figure 1: Connected (a) and disconnected (b) contributions to the scalar density amplitude 
of a nucleon. Please note that all quark lines, including the quark loop, are connected by 
infinitely many gluon lines and virtual quark loops. 

of pion- nucleon scattering data at the unphysical Cheng-Dashen point. A careful analysis 
of the extrapolation procedures has been performed by Gasser, Leutwyler and Sainioflj] by 
means of dispersion relations and chiral perturbation theory. They found a n ^ — 45MeV. 
A consistent value, a n N = 48 ± lOMeV, has been obtained more recently in the framework 
of heavy baryon chiral perturbation theory by the authors of ref. |2j . 

Additional information can be drawn from the baryon octet mass splittings. The flavour 
octet quantity 

a = m " ± md (N\uu + dd- 2ss\N) (2) 
_ 2(N\ss\N) 

anN ~ 1-y ' V ~ (N\uu + dd\N) ■ [6) 

The analysis of the baryon octet mass splittings in first order chiral perturbation theory 
yields o~o — 25MeV. It has been pointed out however ||, that corrections due to terms 
oc m s — 1/2 (m u + mi) 2 may enhance this value on a (20 — 30)% level. Relating these 
findings with the results of the ttN scattering analysis leads to the estimate y ~ 0.2 — 0.4. 
This implies that strange quark loops decrease the value of <To by 20 — 40%. Naively, i.e. 
under the assumption of an approximate SU (3) flavour symmetry of the disconnected parts 
of the nucleon scalar density, one would expect light quark loops to contribute similarly to 
the pion-nucleon a-term. 

ct^tv as well as y have been studied recently in quenched lattice simulations by Fukugita 
et al.ff] and Dong et al.||. Both calculations find a value for cr nN consistent with the 
'experimental' results quoted above, and a rather large ratio 
< N\uu + dd\N > disc I < N\uu + dd\N > co „~ 2/1. 

Furthermore, ref. estimates y = 0.66(15). Given the value of 2 : 1 for the ratio of 
disconnected to connected contributions, this is exactly the value one would expect from 
the assumption of flavour symmetry of the disconnected contributions. Such a finding 
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appears plausible since the quenched QCD vacuum is sensitive neither to quark flavour 
nor mass. 

The authors of ref. H on the other hand find a much smaller value, y = 0.36(3). But 
this is mostly due to their use of phenomenologically inspired ansatze for the extrapolation 
in the quark mass and for the renormalization of disconnected contributions. We will come 
back to this point later. Moreover, their data sample being very limited, one might have 
doubts on the reliability of their error analysis. 

Apart from this, it is by no means obvious that quenched lattice simulations are at 
all suited to yield sound first principle QCD estimates for a n N and y. For the quenched 
approximation neglects internal quark loops in the vacuum field configurations. As sea 
quark loops are essential in the calculation of disconnected contributions, the quenched 
approximation appears to be inconsistent and might introduce a serious systematic bias to 
a^N and y. 

It is therefore of utmost importance to study these quantities in full QCD lattice simula- 
tions. Apart from the issue of the numerical value of a^N one might learn about the physics 
of sea quarks in QCD from the relative weight of disconnected to connected contributions 
and the amount of flavour symmetry breaking in the disconnected sector. 

We emphasize that the size of the y parameter is determined by these quantities: a 
large ratio of disconnected to connected parts and a high degree of symmetry breaking 
impose a low value on y. 

The calculation of disconnected correlations functions in lattice QCD requires a high 
statistics of gauge field configurations, since such correlators receive contributions only 
from vacuum fluctuations. Even with a 'state of the art' statistics of 0(200) configurations 
one needs elaborate analysis techniques to enhance the signal to noise ratio. Previous 
exploratory full QCD simulations H, [7J have therefore not been sensitive enough to resolve 
a possible unquenching effect. 

In this paper we present the results of a full QCD lattice simulation with nf — 2 
quark flavours of (mass degenerate) Wilson fermions at (3 — 5.6 and a lattice volume of 
nl x n t = 16 3 x 32 points. This corresponds to a lattice cutoff a" 1 ~ 2.3GeV and a 
spatial length L ~ 1.4fm. We have generated 200 statistically independent vacuum con- 
figurations at each of our 4 values of the sea quark mass, which correspond to m^/m p = 
0.833(3), 0.809(15), 0.758(11) and 0.686(11). Details of the simulation as well as the anal- 
ysis of the light hadron spectrum can be found in ||. 

The paper is organized as follows. In chapter 2 we will define the lattice correlators 
and explain our methods to extract the connected and the disconnected scalar density 
amplitudes of the nucleon. In chapter 3 the raw results are presented and the quality of 
our signals as they emerge from different analysis methods is discussed. The extrapolation 
to light quarks is performed in chapter 4. Here we also obtain our results for cr^N and y. 
Finally, discussion and conclusions are given in chapter 5. 

2 Analysis Setup 



3 



2.1 Ratio Methods 



In order to optimize the signal to noise ratio and to study the systematics of different 
procedures we have applied several analysis methods both to the connected and to the 
disconnected parts of the scalar density matrix element. 



2.1.1 Global Summation Method 

The standard procedure is to consider the ratio 

R SUM (t) = (4) 
ggOTg, 0) Ey, yo [gg] (y, Vo)N(x, tj) 



Eaf<JVt(0,0)iV(f,*)) 



[w] (ViVo)) 



N denotes an interpolating operator for the nucleon. With the help of the generating 
functional formalism of the path integral, and using the Feynman-Hellmann theorem, one 
obtains 

R SUM (t)=A+(N\qq\N)t, (5) 

provided the nucleon is in its ground state. Naively one would expect this requirement to 
hold at large time distances t. However, due to the summation over all time positions of 
qq in eq.f|, R SUM might even then be contaminated by nucleon excitations. Furthermore, 
the ratio receives contributions from quark loops at distances yo » t. Those do not add 
substantially to the signal, but, especially for the disconnected part, might enhance the 
noise. It is therefore advantageous to undo the summation over y by evaluating the signals 
at definite values yo, in the range <C yo <C t. 

2.1.2 Plateau Density Method 

The local density ratio 

R PLA (t,y Q ) = (6) 

^(N^o^EyM (y,yo)N(x,t)) 



Zg(M(p,0)N(2,t)) 



- N] (y>yo)) 



allows for an isolation of the proton ground state with respect to both t and yo- Its asymp- 
totic time dependence can be evaluated in the transfer matrix formalism, without recourse 
to the Feynman-Hellmann theorem. For ^ y ^ t the ratio becomes independent of t 
and yo, and one obtains[| 

R PLA (t,yo) = (N\qq\N) . (7) 

Thus, the ground state signature of R PLA is a plateau instead of a linear rise. The height 
of the plateau will tell us about the scalar density nucleon matrix element. 

1 This form is valid for an infinitely extended lattice in time. 
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We will demonstrate below that this method works well for the connected parts. How- 
ever, one looses the advantage to gain statistics by summing over y , which becomes a 
crucial issue for the disconnected part. For the calculation of the latter we have therefore 
decided to use a slightly modified technique which still accumulates data with respect to 
I/O) but in a region characterized by ground state dominance of the signals. 



2.1.3 Plateau Accumulation Method 

The plateau accumulation method (PAM) combines the advantages of the global summa- 
tion and of the plateau density techniques. It is defined by 

t-At 

R PAM (t,At Q ,At) = £ R PLA (t,y ), (8) 

yo=At 

with 1 < At, At < t. The asymptotic time dependence is given by 

R PAM (t, At , At) = B+ (N\qq\N)(t -At- At ) . (9) 



2.2 Numerical Evaluation 
2.2.1 Connected Contributions 

We compute the numerator of eq.|] (summation method) with the standard insertion tech- 
nique 0. This method is advantageous if one has to sum over all time positions y . The 
additional effort, on top of the standard quark propagator calculation, is just to compute 
a modified quark propagator A(x, 0), defined as the solution to 

MA(x, 0) = A(x, 0) . (10) 

Here M is the fermion matrix and A the quark propagator. Graphically, the modified 
quark propagator corresponds to the quark line with a cross in fig. |l]a. 

The standard insertion technique can be in principle applied also to the numerator 
of eq.|6] (plateau method). However, as one avoids the summation over y in this case, 
one would have to solve eq.|10| n t times. Instead, we use the advanced insertion technique 
proposed by Martinelli and SachrajdaJTIJ]. Here, one keeps the nucleon sink and source 
at the largest possible time separation^, and computes an advanced propagator by solving 
the equation 



M{ N(0) ^———^ N(x) } = N(0) ^——— ^N(x) ■ ^ 

Note that the r.h.s. is just a combination of standard quark propagators. The 3-point 
correlation is then given by the product of the advanced and the standard quark propagator. 
The position in time, yo, of the quark density qq is not fixed here and can be varied without 
additional cost. 



2 For (anti)periodic boundary conditions, the largest possible time distance is n t /2. 
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2.2.2 Disconnected Contributions 



The determination of the quark loop contributions to the scalar density matrix element 
requires the calculation of the trace of the quark propagator 

L (Vo) = W\% 2/o, a, a; V, 2/o, a, a) = TrA(y ; y ) . (12) 

y,a,a 

a and a denote Dirac and colour degrees of freedom. An exact determination of L with 
conventional Krylov subspace methods would be prohibitively expensive since one would 
have to apply such an algorithm x n t times. Instead one has to rely on approximate 
methods like the volume source technique[TT] and the stochastic estimator technique with 
Gaussian[0 and with ^[l^] noise. The latter allow to control the accuracy of the esti- 
mator on each gauge configuration. 

It has been demonstrated recently |TJj] that the stochastic estimator technique with 
(complex) Z2 noise is superior for our lattice setup. Therefore we used this method here, 
albeit with small modifications which allow for the determination of 

Hvo, ft, ft) = 2 A( #' y°> a ' z 3 '' v> y°> a > P) > ( 13 ) 

y,a 

where the Dirac indices are not contracted. Thus, the results can be re-used to calculate 
vacuum loops other than scalars. For completeness we sketch the relevant formulae of this 
spin explicit method here. 

On a given configuration and for each estimate we choose iV^ complex Z 2 random 
vectors r](y, yo, a, a), with x n t x 4 x 3 entries. Each component of r] has the properties 

r)*(i)r)(i) = 1 , (v*(i)v(j)) = for i ^ j . (14) 

The brackets of the right equation denote the average over (infinitely many) stochastic 
estimates. From 77 we compose 4 spin explicit random vectors 

V P (y,yo,a,at) = r}(y,yo,a,a)5 at/3 ,/3 = 0,1,2,3. (15) 

Note that there is no sum over a on the r.h.s. . The vector An^ is then obtained as the 
solution to 

M[An p \ = if . (16) 

According to eq.[T3, the product 

P(t ,(3',f3) = (f'^A^, (17) 
with ^'' i0 (y, y , a, a) = rf(y, y Q , a, P)5 y0ito 5 a>p , 

converges in the limit N E — > 00 on each gauge configuration to 

(P(t ,(3',(3)) = L(y ,(3',(3). (18) 

In this work we have used 100 stochastic estimates per configuration. We have checked 
that this suffices to reach the asymptotic region. 
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Figure 2: Summation method: The raw data R u and Rd for the connected amplitudes C u ^ 
at our sea quark masses. The fits (range and value) are indicated by solid lines. 



3 Raw Data 

Figs.0 and |3] display the ratios R SUM and R PLA of the connected contributions for our four 
quark masses. On each graph we show the data for both, the (scalar) interactions of the u 
and of the d quarks of the proton. 

Note that R PLA is plotted in figJS clS £L function of the (time) position of the scalar 
interaction y Q . The time separation of proton sink and source is fixed at t = 16. 

All ratios exhibit clear signals, even for the smallest value of our quark masses. From 
fits according to eqs.|5] and [7] we extract the connected parts of the scalar density matrix 
element of the proton. These are listed in tabjl]. 
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Figure 3: Plateau method: The raw data R u and Rd for the connected amplitudes C U) d at 
our sea quark masses. The fits (range and value) are indicated by solid lines. 



Both methods yield consistent results within statistical errors. It appears however, 
that the summation method systematically leads to slightly larger values than the plateau 
method. We attribute this to a small contamination of the SUM results with excited 
proton contributions. We will therefore use the data from the plateau method for our final 
analysis. 

The signals of the disconnected contributions are shown in fig. ^ for the summation 
method and in fig. [5] for PAM. 

It turns out that the data analysed with the summation technique is quite noisy and the 
signal tends to vanish at the smallest quark mass. In contrast, PAM produces significantly 
improved signal to noise ratios over our entire range of quark masses. 
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Method 


K 


r 


Ch 


r 

' u-\-d 




0.1560 


2.08(4) 


1.13(2) 


3.21(5) 


SUM 


0.1565 


2.13(6) 


1.18(4) 


3.31(9) 




0.1570 


2.26(11) 


1.34(6) 


3.60(17) 




0.1575 


2.40(16) 


1.46(13) 


3.86(28) 




0.1560 


1.92(5) 


1.05(4) 


2.98(8) 


PLATEAU 


0.1565 


1.99(8) 


1.15(5) 


3.15(13) 




0.1570 


2.08(13) 


1.27(11) 


3.35(24) 




0.1575 


2.26(14) 


1.33(11) 


3.59(23) 



Table 1: Lattice results (unrenormalised) of the connected amplitude C q = 

\P\ K aea) \QQ\ K sea) |-P( K sea))con- 



The outcome of the fits according to eqs.[5| and |] is compiled in tab.^|. The PAM data 
are clearly superior with respect to statistical errors. It is gratifying to find the results of 
both methods to be compatible within (large) statistical uncertainties. 

In order to assess the systematics of PAM we have evaluated R PAM for several values 
of At = At. As can be seen from fig. || D q shows no dependence on At, within our our 
statistical precision (200 gauge configurations). Throughout this paper we have therefore 
used At = At = 1. 

In tab.0 we have also included the raw data for the strange loop matrix element 
(P\ss\P), which will be used for the determination of the y parameter. To obtain this 
data we held the loop quark fixed at the strange quark mass (c.f. M), and kept valence 
and sea quark masses degenerate under chiral extrapolation. 

4 Physics Results 
4.1 Pion-Nucleon cr-Term 

Our procedure to obtain the physical value of a n N is to extrapolate the results of tabs.|I], 
H with respect to the quark mass to m u d, and to multiply subsequently with the lattice 
cutoff a" 1 and m uc i. Note that no renormalization is necessary, as a^N is a renormalization 
group invariant quantity. 

In fig.[7] we display the extrapolations of the connected and disconnected amplitudes. 
Since the statistical quality of the disconnected contribution does not allow to resolve for 
higher orders in m q , we have decided to use consistently a linear ansatz for all contributions. 
We emphasize that this is equivalent to a quadratic ansatz in M^(m q ), since the Feynman- 
Hellmann theorem yields dM N /dm q (m u d) = (P\uu + dd\P) . 
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The results of the extrapolations are collected in tab|| Note that the ratio of discon- 
nected to connected contributions at the light quark is mass 

(P\uu + dd\P) disc 2D q 

Rd ' c = /pi- ^aa\p\ — = n = L26 57 • ( 19 ) 

(P\uu + dd\P) con C u + C d 

As a cross check we compare the total value C u + Cd + 2D q = 8.82(2.52) with our 
previous result from the quadratic extrapolation of the nucleon mass[§] 

^(rn ud ) = 11.7(4.9) , (20) 
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Figure 5: PAM with Ato = At = 1: T/ie raw date _R 9 for the disconnected amplitudes D q 
at our sea quark masses. The fits (range and value) are indicated by solid lines. 



and find agreement. Note that the statistical uncertainty of the total amplitude obtained 
with the direct (ratio) method is smaller than the one from the derivative method (29% 
compared to 42%). 

Finally, we determine the pion-nucleon a-term in physical units. Witf|]m u d = 0.000901(54) 
(in lattice units) and a^ 1 = 2.30GeV f§ one obtains 

a wN = a - l m ud {P\uu + dd\P) = 18(5)MeV . (21) 

The statistical error has been determined with a full jackknife analysis of the product, 
including the jackknife distributions of m U( i and aZ 1 . 

3 We use the standard definition of the quark mass m q — |(i — -^-). 
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Figure 6: PAM: The dependence of disconnected amplitudes D q on At. 



Clearly, this is a rather small value, compared to the estimate from experiment. The 
latter is quoted to be a n N — 45MeV, which is in rough agreement with quenched lattice 
estimates. 

We emphasize that the origin for this apparent dramatic unquenching effect on a n N 



appears to be the substantial drop observed in the light quark mass m u d- Our result |L6 
m ud = 2.7(2)MeV, which has been obtained using tadpole improved renormalizationfl~5|l 



2k V Ak c 



1 - 0.0098a 



1, 



MS\ 



(22) 



where a A j s (^) = 0.215, is lower than the corresponding quenched estimate at equal cutoff 
(cr 1 ~ 2.3GeV), m u d ^ 5MeV|T9|, [2l| , by roughly a factor of two. 
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Method 


K 




1 y S 


SUM 


0.1560 


0.86(54) 


0.79(54) 




0.1565 


1.76(76) 


1.74(74) 




0.1570 


3.22(1.68) 


3.19(1.61) 




0.1575 


2.09(1.64) 


1.85(1.44) 




0.1560 


1.63(37) 


1.52(35) 


PAM 


0.1565 


1.49(45) 


1.52(44) 




0.1570 


1.67(71) 


1.64(68) 




0.1575 


2.68(93) 


2.49(86) 



Table 2: Lattice results (unrenormalised) of the disconnected amplitudes D q = 
(P(K sea )\qq(K sea )\P(n sea )) dis and D s = (P{K sea )\qq(K str )\P(K sea )) dis . n str = 0.15608. 



Table 3: Lattice results for connected and disconnected amplitudes at Ku g ht = 0.158462. 
We have used the data from the plateau method PAM for the connected and disconnected 
amplitudes respectively. 



C u 


c d 


2D g 


C u + C d + 2D q 


2.40(19) 


1.53(15) 


4.95(2.36) 


8.82(2.52) 



Gupta has reasoned in the context of his discussion on light quark masses |L7] that 
the validity of the one loop formula, eq.|22], is questionable in full QCD, insofar large non 
perturbative contributions due to the presence of sea quarks could arise and readily change 
Z m = Zg 1 by a factor of two, thus compensating the above factor 1/2 in m ud . However, 
this doubt does not provide a source for a possible underestimate to a^-N, since we are 
dealing here with a renormalization group invariant quantity. 

On the other hand, one might question the standard definition and switch to the lattice 
quark mass as defined by use of the axial vector Ward identity |L7|, p0| , [2~l|j . But from the 
final analysis of the CP-PACS group on quenched QCD reported recently one would 
gather a decrease of m ud and a n N by 20 — 50%! One should remember, though, that both 
definitions of m ud differ at our value of a equally from their common continuum limit, and 
it is thus not a priori obvious which one of them is more suitable. 

Obviously the only way to decide on this issue is to perform a scaling analysis of m ud , 
i.e. to repeat its calculation at several lattice cutoffs a -1 , and then to extrapolate to 
the continuum a — > 0. We mention that the CP-PACS collaboration has launched such 
a scaling study in full QCD, based on three different values of the cutoff [113. Their very 
preliminary analysis strongly hints at a small value of the light quark mass in the continuum 
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Figure 7: Linear extrapolations of connected and disconnected amplitudes to the light 
quark mass. Bursts indicate the results of the extrapolation. 

limit, which incidentally is close to the estimate underlying our present work[|16||. 
4.2 y Parameter 

The determination of the y parameter, defined in eq.|3|, requires the calculation of the 
(purely disconnected) matrix element (N\ss\N). In principle, a fully consistent QCD lattice 
analysis of this quantity can be done only in the setting of a realistic rif > 3 simulation, 
in which strange and light quarks enter the dynamics through the fermion determinant. 
Given our nf = 2 setup, we can only resort to a semi-quenched analysis of the strange 
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sector. This is done by identifying the masses of the sea quarks with those of the light 
valence quarks of the nucleon. The mass of the strange quark, which has no counterpart in 
the sea, is contained to match the strange hadron spectrum in semi-quenched analysis^]. 
This procedure accounts at least for the influence of light sea quarks on the above matrix 
element and is certainly more adequate than previous quenched (n/ = 0) calculations. 

The lattice raw results D s = (P(n sea )\ [ss] (n s )\P(K sea )) are listed in tab.||. The values 
differ only slightly from the raw data D q , indicating that an approximate SU(3) flavour 
invariance of the disconnected contributions is still realised in our rif = 2 simulation. 

We extrapolate D s as a function of K sea linearly to the light quark mass. To obtain the 
y parameter we renormalise our result 

2(P\ss\P) = 4.94(2.14) (23) 

by Z s , eq.[22|. This yields 

2Z s (k s )(P\ss\P) 
V = ry } \\ „ - ,-, m = 0.59(13) . (24) 



It has been suggested by Lagae and Liu that the lattice quark mass dependent 
part of the renormalization factor would be quite different for connected and disconnected 
amplitudes. According to their reasoning, based on first order lattice perturbation theory, 
they would lower our value of y by about 30%. We do however not agree with their 
arguments, at least in the case of scalar quark loops, for the following reason: Since the 
combination m q {P\qq\P) is a renormalization group invariant quantity, both contributions 
to the amplitude, connected and disconnected, will be renormalised equally, i.e. by one and 
the same factor Zs = 1/Z m , where Z m renormalises the quark mass. Choosing different 
factors, Z c s on ^ Zf sc , as they propose, would be in conflict with the renormalisation group 
invariance. 



5 Discussion and Conclusions 

We have studied connected and disconnected contributions to the proton scalar density 
amplitude in a nf = 2 full QCD simulation at fixed lattice cutoff and volume. 

It turns out that conventional ratio methods (summation and plateau) yield excellent 
signals for the connected parts, but fail in the determination of disconnected contributions. 
The latter situation can be improved substantially by use of PAM. It is likely that stochastic 
estimator techniques will perform even better on larger lattices, where we would expect 
self averaging effects to reduce the errors from the present level, which is 30 - 40%. 

We find the ratio of disconnected to connected contributions, Rd/ C = 1.26(57), at m ud to 
be lower than the estimate from the quenched simulation at (3 — 5.7 of ||], where they quote 
Rd/c = 2.23(52). Albeit (3 = 5.7 corresponds to rather strong coupling, this difference in 
Rd/c might indicate an unquenching effect. In order to consolidate this tentative conclusion, 
one should (a) repeat the analysis on a larger lattice, where self averaging will help to reduce 
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the statistical errors, and (b) carry out a quenched reference simulation at equal lattice 
spacing and statistics. Work along this line is in progress. 

Due to the small value of m u d in full QCD, the physical result for the pion-nucleon a- 
term comes out rather low. It remains to be seen whether this is a finite cutoff effect. In this 



context it is very interesting that the preliminary results of the CP-PACS collaboration] 19 



obtained in full QCD with an improved action, point at an increase of m U( i by about 30% 
when changing from the standard to the Ward identity definition, at our value of a. It is 
obvious that a scaling analysis of the quark mass in full QCD is of utmost importance. 

So far the y parameter has only been determined with a semi-quenched treatment of 
light quarks. The result, y = 0.59(13), is far away from the phenomenological expectation, 
y ~ 0.2 — 0.4. A full QCD calculation of y with rif > 3 dynamical flavours is therefore 
highly desirable. 
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